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We propose parametrizations of the neutrino mixing matrix. We discuss several Ansatzes: in 
some of them V e 3 = but in others V e 3 typically "comes out" to be of order ~ 0.05. 

I. Neutrino Mixing Matrix 

Wolfcnstein's parametrization [1] of the quark mixing matrix, while not based on any deep theoretical foundation, 
~£f*)~ • has proved to be enormously useful in making sense of the data. Something similar is needed to make sense of 
the neutrino mixing matrix. In reporting the value of a typical mixing angle 9, experimentalists quote variously 
sin#, sin 29, tan#, or their squares. It would be helpful to have a more systematic way of quoting experimental 
values, analogous to what the Wolfenstein parametrization has provided for the quark sector. In the literature, two 
i— H \ parametrizations [2,3] of neutrino mixing had been proposed that I know of. Here we would to propose a quite 
different parametrization. We also suggest a more systematic way for experimentalists to quote their data. As in 
the Wolfenstein parametrization, our discussion is not based on any deep theoretical foundation, but rather on a 
phenomenological analysis of the available data. This paper is thus a somewhat old fashioned mix of Ansatz making 
' and data fitting. 

Thanks to heroic experimental efforts, the neutrino mixing matrix is now known to take on the value 

> : 

CO . / 0.72 - 0.88 0.46- 0.68 < 0.22 

|^|= 0.25 - 0.65 0.27- 0.73 0.55 - 0.84 | (1) 
\ 0.10 - 0.57 0.41 - 0.80 0.52 - 0.83 

, I am citing the report [4] of Gonzalez-Garcia, based in part on the analysis given by Bachall, Gonzalez-Garcia, and 
Pena-Garay [5]. Here the notation |V| is such that the ij matrix element of the matrix \ V\ is equal to \Vij\. 

The mixing matrix V relates the neutrino current eigenstates (denoted by v a (a = e, fx, r) and coupled by the W 
bosons to the corresponding charged leptons) to the neutrino mass eigenstates (denoted by (i = 1, 2, 3)) according 
I ' to 

^ i /'„ 1 I i r, 1 (2) 

We will take the neutrinos to be Majorana [6] as seems likely, so that we have in the Lagrangian the mass term 

C = —v a M a pCvp + h.c. (3) 

where C denotes the charge conjugation matrix. Thus, the neutrino mass matrix M is symmetric. 

For the sake of simplicity we will assume CP conservation so that M is real. With this simplification, M is 
diagonalized by an orthogonal transformation 

/mi \ 

V T MV = m 2 (4) 
V m 3 J 

It will become clear to the reader that everything in this paper is easily generalized to the case in which CP is violated 
so that V is unitary rather than orthogonal, as briefly indicated in section V. Also, we are free to multiply V on the 
right by some diagonal matrix whose diagonal entries are equal to ±1. This merely multiplies each of the columns in 
V by an arbitrary sign. Various possible phases have been discussed in detail in the literature [7,8]. 

II. An Ansatz for the Neutrino Mixing Matrix 

We could suppose either that the entries in V represent a bunch of meaningless numbers (for example, possibly 
varying from domain to domain in the universe) or that they point to some deeper structure or symmetry. In the 
latter spirit, let us make a guess of what V might be. 
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Since V e s appears to be small, let us boldly set it to 0. Next, since l/y/2 
(0.55 - 0.84) = 1 /v/2. Finally, since 1/ V3 - 0.577 we will set V e2 = (0.46 - 0.68) = 
that we know the upper triangular entries of the matrix V : 
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' 0.707 we will guess that — 
l/y/3. In other words, we propose 



(5) 



where X denotes an unknown quantity. 

Remarkably, this essentially fixes the mixing matrix V. Once we take the last column to be proportional to (0, 1, — 1), 
orthogonality and our "knowledge" that V e 2 is 1/V3 immediately fix the second column to be proportional (1, 1, 1) 
and hence the first column to be proportional to (—2, 1, 1). X. G. He and I therefore arrived at the Ansatz or guess 
[9] 
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(6) 



This Ansatz was proposed earlier by Harrison, Perkins and Scott [10]. Also, this mixing matrix (but curiously, with 
the first and second column interchanged) was first suggested by Wolfenstein more than 20 years ago [11] based on 
some considerations involving the permutation group S3. It has subsequently been studied extensively by Harrison, 
Perkins and Scott [12], and by Xing [13]. Attempts to derive this mixing matrix have been discussed by Low and 
Volkas [14]. We will refer to the matrix in (6) as V unless confusion could arise in which case we will refer to it as 
Vhpshz. ■ 

The signs in V were chosen to minimize the number of minus signs. The vanishing of V e s indicates, if this Ansatz 
is correct, that there is no observable CP violating effects in neutrino oscillation experiments. See (26) below. 
The mixing matrix V may be factorized as 



V = VnV, 



23 v 12 



where 
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and 
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(9) 



Note that our sign choices in (6) is such that det V23 = det V\2 — —1 (but det V = 1.) 

The rotation of the neutrino current eigenstates appear to "occur in two steps." In other words, if we follow 
Wolfenstein and define v x = {v^ + v T )/\/2 and v y = (y^ — v T )/^/2 we find that the mass eigenstates are given by 



(10) 



and 



V2= Vf" e + V3^' 



= Vy 



(11) 



(12) 



The matrix V12 describes a rotation through arcsin(-^) ~ 35°. 
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Giunti [15] has proposed another Ansatz, to which we will turn later. 



III. Parametrizing the Neutrino Mixing Matrix 

In order to compare theory with experiment, we have to understand what the experimental values in (1) mean. The 
quoted range of numbers are extracted from different experiments, but the error ranges are obviously not uncorrelated 
since V has to be unitary, or orthogonal with our simplifying assumptions. 

There is not a unique way to extract an orthogonal matrix V cxp from (1). Of the many possible ways we choose 
the following procedure. We take the mean value of the cited range of the various entries in |V| (e.g. for \V e 2\ we 
set (0.46 — 0.68) — ► 0.57 ), fix the signs of various entries to agree with the choices in (6), and thus arrive at the 
"experimental mean values" 

0.00 

Kncan = I 0.45 0.5 0.695 1 (13) 




-0.675 

The matrix Vmean is not orthogonal. Since I^mean^mean is real symmetric it can be diagonalized by an orthogonal 
matrix S so that the matrix K 2 = S T V^ ean V me anS is diagonal with positive definite diagonal elements. We then 
define the "corrected experimental mixing matrix" to be the orthogonal matrix 



V™* = K 1 ean(^- 1 5 T ). (14) 

Numerically, this comes out to be 

/ -0.814 0.578 0.056 \ 

Kxp" = °- 436 °- 545 °- 716 (15) 
\ 0.384 0.607 -0.696 / 

Note that a non-zero value of V e 3 is "generated" but one order of magnitude smaller than the other entries. Second, 
note that V™" is well within the range quoted in (1). Finally, compare V™" with the theoretical Ansatz (6) 

/ -0.817 0.577 0.000 \ 
V= 0.408 0.577 0.707 (16) 
\ 0.408 0.577 -0.707 / 

Clearly, V and V™* 1 are rather close. Thus, three obvious ways of parametrizing the data suggest themselves. We 
could parametrize the data by a matrix W defined by V™* 1 = VW, or a matrix X defined by V™* 1 = XV, or a matrix 
Y defined by V^° p — Vi%YV\i (as suggested by (7)). Namely, we parametrize the deviation of V from V^°" . 

In each of these three possible ways, the matrix (W, or X, or Y) will be close to the identity matrix, and we could 
then follow Wolfenstein and parametrize them accordingly. In each of these cases, we will write the orthogonal matrix 
W, or X, or Y in the form e A ~ I + A + \A 2 with 




A= | n > ) (17) 
We will discuss each of these possibilities in turn. 



(I) We solve the 3 equations Wjj — {V T V e ™ TT )jj, j — 1, 2, 3 (no repeated summation convention) for the 3 unknowns 
Oi,(3, and 7. To the order indicated, we have 3 coupled quadratic equations and we obtain 8 solutions corresponding 
to various choices of signs. We then compute y^yj (Wjj — (V T V,?°")ij) 2 and determine the solution which minimizes 
this quantity. In this way, we find that the best fit to W is given by a = -0.00249, (3 = 0.0440, and 7 = 0.0374. 
The main outcome of this simple phenomenological analysis is that a is an order of magnitude smaller than (3 and 
7, which appear to be comparable. Indeed, we see that a is close to — 3/3 2 /2. In the spirit of Wolfenstein's paper, 
without any theoretical understanding whatsoever of the values of a, (3, and 7, we could try, instead of a 3-parameter 
fit, a 2-parameter fit by setting 

/0 -3v[3 2 /2 -0\ 

A= Zvf3 2 /2 (3 (18) 

\(3 J 
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where we expect v to be of order 1. Solving the 2 equations W\ 2 — (V T V < ?°")i2 and W23 = {V T V™")23, we nn d that 
v = 0.908 and (i = 0.0439. 

(II) The analysis clearly proceeds in the same way: we solve the 3 equations Xjj = {V™"V T )jj and obtain the 
best fit solution a = 0.0413, (3= -0.0143, and 7 = 0.0378. 

(III) In this case we solve the 3 equations Yjj = {V 2 %V™p T V£)jj and obtain the best fit solution a = 0.00249, 
(3 = 0.0143, and 7 = —0.0559. In other words, we have to rotate v ei v x = (17, + v T )j\[2 and v y = iv^ — v T )j\/2 by a 
tiny amount before mixing them with V^. 

IV. Another Ansatz for the Neutrino Mixing Matrix 

Giunti [15] has proposed another interesting Ansatz for the mixing matrix 
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Vg - I 2^2 172 7=2 I (19) 
which also factorizes nicely Vq — V23VG12 where 

V G12 = I 1 4 () i (20) 




Thus, Vg and V differ by a rotation of about 5° in the (1 — 2) plane: 



J_ + 1 I L 

v/2 + 2V3 2 y'e u 

1 



Numerically, 



-0.866 0.5 
V G = \ 0.354 0.612 0.707 ) , (22) 
0.354 0.612 -0.707 

which is to be compared to V in (16). We see that Vg makes the /zl and rl entries smaller (0.354 rather than 0.408) 
and the /i2 and r2 entries larger (0.612 rather than 0.577). 

We could ask whether V or Vq is closer to the experimental mixing matrix V^°p r - The group invariant distance 
between two group elements 0\ and 2 of 5*0(3) is given by 

dia = 3 - Tr(Oj0 2 ) (23) 

We find 

d(V) = 3 - Tr{V T V^) ~ 0.00334 (24) 

and 

d(V G ) = 3 - Tr(V^K c X p r ) ^ 0.0122, (25) 

about a factor of 4 larger than d{V). If (and only if) we believe in the more or less arbitrary way in which V^°" is 
extracted from the raw data (1), then V is a somewhat better fit to the data than V G . 



V. Angular Parametrization of the Neutrino Mixing Matrix 



We could of course always parametrize the neutrino mixing matrix with three Euler angles and a phase. The 
factorization in (7) suggests that we use the analog of the Chau-Keung [16] parametrization commonly used for the 
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Cabibbo-Kobayashi-Maskawa matrix in the quark sector and thus we write (where as usual c 2 3 = cos0 2 3, s 23 = sin 6*23, 
and so forth) 



Kmgular = V 2 ^V Zl V 12 (26) 

10 \ /c 3 i s 3 ie-^\ / -Ci 2 S12 l\ 

C 23 S 23 1 S12 C12 (27) 

s 23 -c 23 / V -s 3 ie i0 c 3 i / V 1 / 
-C31C12 C31S12 s 3 ie - ^ \ 

S12C23 + ci2S 2 3S3ie^ ci 2 c 23 + si2S23S3ie^ S23C31 (28) 

S12S23 - Ci2C 2 3S3ie^ Ci 2 S 23 + Si 2 C 2 3S 3 ie^ -c 23 c 3 i J 

(Our parametrization differs slightly from that of Chau and Keung in that we take det V23 = det V\ 2 = — 1 in 
accordance with our earlier sign choice.) Evidently, for the Ansatz in (6) and in (19) 9 3 i — and thus V e s = 
ssie - ** = 0. 

If we take V^£ r seriously, then a fit to V an g U iar gives 

012 = 0.617 = 35.4°; 23 = 0.800 = 45.8°; 3 i = 0.056 = 3.21°. (29) 
These values are extremely close by construction to the corresponding values for the Ansatz Vhpshz in (6) 

012 =35.3°; 23 = 45°; 03i = 0°. (30) 



This discussion makes clear how easy it is to construct more theoretical Ansatzes. Denote the three columns in the 
theoretical V by v\ , v 2 , and V3 respectively. Pick one of the columns, call it v, compare with the corresponding column 
in V™p T and fit two of the components in v to some combinations of (due to some possibly unjustifiable theoretical 
prejudice) small integers and their square roots. The third component of v is then fixed by orthonormality. Pick 
another column, call it if, and fit one of the components in v' to some combinations of small integers and their square 
roots. The other two components of if are then fixed by orthonormality and by orthogonality to v; they emerge as 
solutions of a quadratic equation. The last column is then fixed to be / = ijx if up to a sign. 

V. To Make Some Experimentalists Happy 

A number of major experiments to detect CP violation in neutrino oscillation are being planned. Unfortunately, if 

V e 3 vanishes the observable signal for CP violation also vanishes, as is shown clearly in (26). The vanishing of V e z, 

however, is the one common feature of Vhpshz and Vg- Thus, the experimentalists involved have a human tendency 
to strongly dislike both Vhpshz and Vg- 

Suppose, in order to make some experimentalists happy, we try to construct an alternative Ansatz Vh with (Vh )e3 = 
1/5, close to its experimental upper bound (1) of 0.22. Of course at each step of the construction of Vh we will have 
to make an arbitrary choice. We will choose what I consider to be the most reasonable. Suppose we want (Vj?)^ 
and (Vh)t3 to remain equal in magnitude. This fixes the last column of (Vh) to be {1, 2-^/3,-2-^/31/5. Note that 
(Vf/) M 3 = 2\/3/5 remains very close to 1/V2, its value in both V and Vg, since its square is equal to 12/25. Let 
us next suppose (Vh) E 2 = l/\/3 as suggested by experiment. Orthogonality then fixes (Vh)^ and (Vh) T 2 to be 
(V141 — v / 3)/12v / 3 and (V141 + v / 3)/12v / 3 respectively. The first column is completely fixed to be the vector cross 
product of the second and third columns. We obtain what some might consider a rather "ugly" matrix 



V H = 



( _I /iz J_ I \ 

5 V 3 V3 5 1 

5 I V47 VT41-V3 2V3 

12 60 12^3 5 

\ 12 60 12^3 5 / 



(31) 



Theoretically, even if we were optimistic enough to imagine that we could derive V or Vq (from some discrete symmetry 
group [14] for example) it is hard to imagine how we could possibly derive Vh- However, Vh provides a reasonable fit 
to experiment by construction; numerically 



5 



-0.792 0.577 0.200 
V H = ( 0.531 0.488 0.693 
0.302 0.655 -0.693 



(32) 



to be compared with (15). We note that 



d(V H ) = 3 - Tr(VgV£?) c 0.0214 



(33) 



about a factor of 2 larger than d(Vc)- This somewhat meaningless exercise merely shows that it is perfectly consistent 
with present data to have a large e3 matrix element, but it also illustrates the fact, which the reader can convince 
himself or herself by playing around, that it is not easy to come up with a "nice" theoretical Ansatz for the mixing 
matrix. Orthogonality typically generates square roots of large integers, and by going through a number of cases, I 
see that this usually happens unless V e3 = as in Vhpshz, as may be obvious to the reader. 

VI. The Value of V e3 

As mentioned earlier, the value of V e3 is of crucial interest to some experimentalists. In the last section we simply 
put in the value 1/5, and in both the theoretical Ansatzes (6) and (19) V e3 is set to 0. We could perhaps turn things 
around and try to "predict" V e3 . 

We will try to find another Ansatz, which we will call Vc, by fitting 3 numbers from the first and second columns 
of V^°" , and see what we get for V e 3- For the reader's convenience I reproduce V^°" here: 



T/corr 
*exp 



-0.814 

0.436 

0.384 



0.578 0.056 
0.545 0.716 
0.607 -0.696 



(34) 



Let us look at column 2. As I already noted in constructing (6), ~ 0.577 is very close to (V™") e 2 and so we set 
{Vc)e2 = -yg- But in contrast to the construction of (6), here we will take seriously the apparent fact that (V™ p ")n2 

is somewhat smaller than (V r ™ p rr )e2 : with a bit of fiddling we find that §^/§ — 0.544 gives a good fit to (V^p 1 )^- 

Orthonormality then fixes (V^°") T 2 to be — 0.609. Now that we have filled column 2 of Vc we can only fix 

one of the numbers in column 1: orthonormality of column 1 and orthogonality between column 1 and column 2 
will fix the other two. We observe that \\J\ — 0.385 is very close to {V^° p TI ) T i, so let's set (Vc) T i = We 

now have to solve a quadratic equation to determine (Vc) e i and (Vc)^i, and here comes a "problem": as in the last 
section, there is no reason that the solution of a generic quadratic equation would involve only small integers and the 
square root of small integers. In fact, already \/T0 has appeared. Indeed, we obtain (Vc)ei = — -^(3^26 + V30) and 
{Vc)ni — I53(27\/13 — 8V1~5). The vector v 3 in column 3 is now determined to be the vector cross product of the 



vectors v\ and V2 in column 1 and column 2, and we obtain the Ansatz 



/-A (3 V26 + V30) ^ #(-12 + 



^(27^13-8^) 



i 

V3 

2 PI 

3 V 3 

_L 

3 



T§3( 27 + 2^ 



/l95)\ 
/ T95) 



(35) 



Thus, we "predict" V e3 to be 



V2, 



— (V195 - 12) ~ 0.0544677, 
51 

which "happens" to be in a range encouraging to experimentalists. Numerically, 

-0.815 0.577 0.055 
V c = ( 0.434 0.544 0.718 
0.385 0.608 -0.694 



(36) 



(37) 



By inspection we see that (37) gives a good fit to (34 ); this is confirmed by the group invariant distance between Vc 



and K™": 
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d(V c ) = 3 - Tr(VgV c c x ° p rr ) ~ 1.08 x 10~ 5 . 



(38) 



Compare d(Vc) with rf(l^) and d(Vc)- 

Of the three Ansatzes (I am discounting Vh) Vc emerges as the champion. However, most theorists, perhaps 
unjustifiably prejudiced against large integers (they might pop out of string theory?), would probably regard Vc as 
less elegant than the other two Ansatzes even though it gives a better fit. Thus, in a sense, this exercise reinforces 
my belief in, or rather, liking for the Ansatz V in (6). 

Whenever experimentalists refine the raw data in (1) the reader can easily play this game. Extract some kind of 
mean value V mcan and then orthogonalizc or unitarize it to an orthogonal or unitary V^?p r - The proposed Ansatz 
should be judged by its group theoretic distance to V^°" . This procedure is further illustrated in the Appendix. 



VII. Neutrino Mass Matrix 



Neutrino oscillation experiments can only determine the absolute value of the mass squared differences Am?- = 
mf — m 2 rather than the three individual masses m,. At the 99.3% confidence level Am?- are determined to be 
1.5 x 10" 3 eV 2 < |Am§ 2 | < 5.0 x IQ-^eV 2 and 2.2 x 10" 5 eV 2 < |Am^| < 2.0 x lQ- A eV 2 , with the best fit values 
given by 

|Am 2 2 | = 3 x l(T 3 eV 2 (39) 

and 

| Am 2 ! | = 7 x 10~ 5 eV 2 (40) 

We could have either the so-called normal hierarchy in which |m 3 | > |m 2 | ~ |mi| or the inverted hierarchy \ms\ < 
|m 2 | ~ \m\\. At present, we have no understanding of the neutrino masses just as we have no understanding of the 
charged lepton and quark masses. 



Suppose wc claim that we know the neutrino mixing matrix V. As before, call the three column vectors in the 
mixing matrix Vf. The neutrino mass matrix is then given by 

3 

M = Y,m l v l {v l ) T . (41) 

i=l 

In other words, the 6 parameters of a real symmetric matrix have been reduced to the 3 eigenvalues {mi}. In particular, 
if we believe in Vrpshz we have 

( A -2 -2\ / 1 1 l\ /0 

In [9], in discussing the neutrino mass matrix, we proposed a basis of 3 matrices other than those that appear in 
(42). First, the 3 column vectors in V are the eigenvectors of the matrix 

2 

M = a \ -1 3 | (43) 
3 -1 

with eigenvalues mi = m 2 = 2a, and m 3 = —4a. (The parameter a merely sets the overall scale.) Thus, with M as 
the mass matrix Am 21 = and this pattern reproduces the data | Am|i | / 1 Am§ 2 1 -C 1 to first approximation. 

Because of the degeneracy in the eigenvalue spectrum, V is not uniquely determined. To determine V, and at the 
same time to split the degeneracy between m\ and m 2 , we perturb M to M = M + SMp, with the " democratic" 
form 

111 

5M D = ea | 1 1 1 ) . (44) 
1 1 1 
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The matrix 5Md is evidently a projection matrix that projects the first and third columns in V to zero. Thus, the 
eigenvalues are given by mi = 2a, ni2 — 2a(l + 3e/2), and 7713 = —4a, where to the lowest order e = Ara^/Am^ and 
a 2 = At7J3 2 /12. Finally, to break the relation \ms\ = 2\mi\ ~ 2|m2| we can always add to M a term proportional to 
the identity matrix. 

VIII. Conclusion 

In this note, we made the simple point that experimental data should be summarized and parametrized by its 
deviation from some suitable theoretical Ansatz, so that we arc effectively dealing with an orthogonal or unitary 
matrix close to the identity, which we can then study in the same spirit as the Wolfenstein parametrization. Here for 
simplicity we have focussed on the orthogonal case, but it is straightforward to generalize our discussion to the unitary 
case by adding a phase. We also advocate using the group invariant distance between two matrices to measure the 
goodness of fit. 
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Appendix 

In the Ansatz Vc discussed in section VI, I first set column 2 in Vc to 

( 4\ 

V2 




2 

3' 
1 

3 ' 



(45) 



/10 J 
3 / 



and V T i to j\J\- I then determined V e \ and V^\ in 



vi 




(46) 



by solving the quadratic equation resulting from the conditions vi ■ v\ = 1 and v\ ■ V2 = 0. Finally, 



«3 




(47) 



was determined by v 3 = v\ x v 2 . "Large" integers and a non-zero value of V e z emerge. 

Now suppose, instead of fixing V T \ \nv\, I follow the two alternate possibilities: (A) fix V e \ to 



-0.817, or 



(B) fix to 



0.433. Note for comparison that (V™ rr ) e i = -0.814 and (V c 



corr\ 
exp J/il 



0.436. For possibility (A) 



and possibility (B), let me ask you three questions. (1) Does a non-zero value of V e s emerge? (2) Do "large" integers 
appear? (3) How close is the resulting Ansatz to V^°"? 

The point is that, at least for mc, it is not easy to answer these questions without doing the computation (it is of 
course easy to write a short program to perform the computation) . 
It turns out that 



/ 



V A 



2 J_ 

3 V3 
2 l~2 

3 3 V 3 

5 1 

3 3" 





VE 
3 



(48) 



- / 

3 / 



has a "nice" form without "large" integers and predicts V e 3 — 0, with a distance from V^°" of 



d(V A ) = 3- Tr{VjV™ T ) ~ 0.00484, 



(49) 
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better than Vc and slightly worse than Vrpshz- Meanwhile, Vb contains huge integers and is too "ugly" to write 
down here, but predicts 

V e3 = -i- (27\/l0 - 2V1338) ~ 0.0536 (50) 

228 

in a range encouraging to experimentalists, and deviates from V™* 1 by only 

d(V B ) = 3 - Tr(V£ Kxp r ) - 1-78 x 10" 5 , (51) 

much better than both Vc and Vhpshz and only somewhat worse than Vc- 

It appears that the absence of large integers and V e3 = are correlated. If we take V^°" seriously and if we arc not 
prejudiced against large integers we can get a better fit with a non- vanishing V e3 . 
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